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1 Introduction 

1.1 Integrable evolutions of space curves 

Much of the work on integrable curve evolution equations has been guided by the fundamental 
role played by the differential invariants of the curve (e.g., curvature and torsion in the Euclidean 
setting) in helping identify the curve evolution as an integrable one. Perhaps the most important 
example in the case of space curves is that of the Localized Induction Equation (LIE) 



describing the evolution of a curve with position vector l{x,t) in M'^, and Euclidean arclength 
parameter x. The complete integrability of equation (1.1) was uncovered by the realization, due 



one of the two best-known integrable nonlinear wave equations (the other being the KdV equa- 
tion) . 

In this paper, we also use as a guiding principle the observation that many (but not all) inte- 
grable curve evolutions have the property of local preservation of arclength, i.e., the associated 



In honor of Peter Giver. 



7t = 7^ X 1., 



(1.1) 




*This paper is a contribution to the Special Issue "Symmetries of Differential Equations: Frames, Invariants 
and Applications". The full collection is available at http://www.emis.de/journals/SIGMA/SDE2012.html 
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vector fields satisfy a non-stretching condition. For example, the LIE vector field W = 1 x ^ ^xx 
satisfies the condition 5vk||7x|| = 0, where 5w denotes the variation in the direction of W . Thus 
the local arclength parameter x is independent of t, and the compatibility conditions Ixt = ^txi 
"^xxt = '^txx, "^xxxt = "^txxx (more commonly written as compatibility conditions of the Frenet 
equations and the evolution equations for the Frenet frame) turn out to be equivalent to the 
Lax pair of the NLS equation for ip. 

Indeed, many integrable curve evolutions in various geometries have been found by looking 
for non-stretching vector fields that produce compatible equations for the moving frame of the 
evolving curve; in the case of space curves, the geometries explored include Euclidean [11], 
spherical [6], Minkowski [21], affine and centroaffine [4]. (Moreover, integrable curve evolutions 
without preservation of arclength have been found in projective [17], conformal [16] and other 
parabolic geometries.) The approach in these investigations involves finding suitable choices for 
the coefficients of the non-stretching vector fields (relative to a Frenet-type frame) and often 
assuming special relations among the differential invariants; thus it can be challenging to identify 
integrable hierarchies. 

Another approach to investigating the relation between a non-stretching curve evolution and 
the integrable PDE system satisfied by the differential invariants is to seek a natural Hamil- 
tonian setting for the curve flow. The LIE was shown by Marsden and Weinstein [18] to be 
a Hamiltonian flow on a suitable phase space endowed with a symplectic form of hydrodynamic 
origin (see also [1,2]). In a fundamental paper [14] Langer and Perline used this framework to 
explore in depth the correspondence between the LIE and NLS equations and, along the way, de- 
rived a geometric recursion operator at the curve level that made it easy to obtain the integrable 
hierarchies of both curve and curvature flows, as well as meaningful reductions thereof [12, 13]. 

In this article we study integrable evolution equations for closed curves in centroaffine 
beginning, as in [14], with a natural pre-symplectic form on an appropriate infinite-dimensional 
phase space. The Hamiltonian setting allows us to construct integrable hierarchies of curve flows 
and the associated families of integrable evolution equations for the centroaffine differential 
invariants (which turn out to be equivalent to the Boussinesq hierarchies). The motivation for 
addressing the centroaffine case comes from an interesting article by Pinkall [26], who derived 
a Hamiltonian evolution equation on the space of closed nondegenerate curves in the centroaffine 
plane. The simple definition of the symplectic form in the planar case (related to the SL(2)- 
invariant area form) suggests that an analogous description may be possible in the 3-dimensional 
case, where a parallel could be drawn with the more familiar Euclidean case treated by [14]. 

Before describing the organization of the paper, we briefly discuss Pinkall's original setting 
and some results of ours for the planar case. 

1.2 Pinkall's flow in 

Centroaffine differential geometry in refers to the study of submanifolds and their proper- 
ties that are invariant under the action of SL(n), not including translations^. For example, 
a parametrized curve 7 : / — )• M" (where I is an interval on the real line) is nondegenerate if 



for all X £ I, and this property is clearly invariant under the action of SL(n). Thus, for these 
curves the integral 



det (7(x),7'(x),...,7("-i)(x)) /O 




(1.2) 



^Some authors [23] refer to this geometry as centro-equi-affine due to the choice of the unimodular group 
SL(n), while using centro-affine to refer to geometry invariant under the general linear group GL(n). 
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is SL(n)-invariant, and represents the centroaffine arclength, where for the sake of convenience 
we use the notation 

l-'^ij ■ • • , Xn\ := det{Xi, . . . , X„) 

for n-tuples of vectors Xi G M". (The fractional power in (1.2) is necessary to make the integral 
invariant under reparametrization.) 

In the case where n = 2, Pinkall [26] defined a geometrically natural flow for nondegen- 
erate curves in M?, which he referred to as star-shaped curves, as follows. Suppose that 7 
is parametrized by centroaffine arclength s, so that |7,7'| = 1 identically. It follows that 
Iss = — p(s)7 where p{s) is defined as the centroaffine curvature. Along a closed curve 7, one 
defines the skew-symmetric form 

ujiX,Y) = I |X,y|ds, (1.3) 

where X and Y are vector fields along 7. This pairing is nondegenerate on the space of vector 
fields that locally preserve arclength. Then the symplectic dual with respect to (1.3) of the 
functional §^p{s) ds is the vector field 

X = \p,l-pl,. 

Pinkall's flow 7t = \ps'^ — P'^s induces an evolution equation for curvature that coincides 
with the KdV equation, up to rescaling. In an earlier paper [3] , we showed how to use solutions 
of the (scalar) Lax pair for KdV to generate solutions of Pinkall's flow. In particular, we showed 
that varying the spectral parameter in the Lax pair for a fixed KdV potential q corresponds to 
constructing a solution to the flow with curvature given by a Galileian KdV symmetry applied 
to q. We also derived conditions under which periodic KdV solutions corresponded to smoothly 
closed loops (for appropriate values of the spectral parameter) and illustrated this using finite- 
gap KdV solutions. 

1.3 Organization of the paper 

In Section 2 we introduce basic notions concerning the differential geometry of nondegenerate 
curves in centroaffine M^, including centroaffine arclength, differential invariants, and non- 
stretching curve variations. This section also contains a discussion of the relation between 
nondegenerate curves and parametrized maps into MP^. In Section 3 we generalize Pinkall's 
setting to by introducing a pre-symplectic form on the space of closed unparametrized star- 
like curves; we also compute Hamiltonian vector fields associated with the total length and 
total curvature functionals. Flow by these vector fields induces evolution equations for the 
differential invariants; we discuss these equations in Section 4, including their bi-Hamiltonian 
formulation. Lax representation, and the connection with the Boussinesq equation. In Section 5 
we show that the Poisson operators introduced in Section 4 give rise to the Boussinesq recursion 
operator, generating a (double) hierarchy of commuting evolution equations for the differential 
invariants. In Theorem 5.4, we relate the Hamiltonian structure for starlike curves and the 
Poisson structure for the differential invariants, and obtain a double hierarchy of centroaffine 
geometric evolution equations. We conclude Section 5, and the paper, by considering which of 
these flows preserve the property that 7 corresponds to a conic under the usual projectiviza- 
tion map tt : — )• MP^. We show that the sub-hierarchy of conicity- preserving curve evolutions 
induces the Kaup-Kuperschmidt hierarchy at the curvature level. 
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2 Centroaffine curve flows in 
2.1 Centroaffine invariants 

Let 7 : / — 7- M'^ be nondegenerate. We parameterize 7 by centroaffine arclength, so that 

|7,7',7"| = 1. (2.1) 

We assume for the rest of this subsection that x is an arclength parameter. 
It follows by differentiating (2.1) with respect to x that 

7"'=po7 + pi7' (2.2) 

for some functions po{x) and pi{x). As explained below, these constitute a complete set of 
differential invariants for nondegenerate curves. 

Remark 2.1. Huang and Singer [10] refer to nondegenerate curves in centroaffine M'^ as starlike. 
They define invariants k and r which correspond to —pi and po respectively. Labeling pQ as 
torsion is appropriate, since nondegenerate curves that lie in a plane in (not containing the 
origin) are exactly those for which pQ is identically zero. 

Remark 2.2. Some insight into the meaning of the centroaffine curve invariants can be gained 
by considering the relationship between 7 and the corresponding parametrized curve T = vr o 7 
in MP^, where vr : M'^ — MP^ is projectivization. The nondegeneracy condition on 7 corresponds 
to T being regular and free of inflection points. Conversely, any such parametrized curve 
T : M MP2 has a unique lift to 7 : M — ?■ which is centroaffine arclength-parametrized; we 
refer to 7 as the canonical lift of T. When written in terms of T instead of 7, the invariants po 
and pi are (up to sign) the well-known Wilczynski invariants [30]. Since these invariants define 
a differential equation whose solution determines the curve uniquely up to the action of the 
group SL(3), any other differential invariant must be functionally dependent on po, pi and their 
x-derivatives. 

According to Ovsienko and Tabachnikov [25], the cubic differential {po — ^p\){dx)^ has the 
interesting property that it is invariant under reparametrizations of T. Curves in MP^ for which 
this differential vanishes identically are conies. For curves for which the coefficient po — ^p'^ is 
nowhere vanishing, one can define the projective arclength differential {po — ^p\)^^^dx. Those 
parametrized curves in MP^ for which pQ — ^p'^ = C (a nonzero constant) are parametrized 
proportional to projective arclength, and we use the same terminology for their canonical lifts 
into MP. (Note that, in this case, the projective arclength differential is C^/^ times the cen- 
troaffine arclength differential dx.) 

Along a nondegenerate curve, an analogue of the Frenet frame is provided by vectors 7, 7', 7". 
In fact, if we combine them as columns in an SL(3)-valued matrix W = (7, 7', 7"), then the 
analogue of the Frenet equations is 

/O po' 

w^ = wii Pi 

\0 1 0, 

However, for later use it will be convenient to define a different SL(3)-valued frame F{x) 
(7, 7', 7" — pil) which satisfies the Frenet-type equation 



F, = FK, K= \1 , (2.3) 




where ki = pi, k2 = Po — p'\- Of course, fci, k2 also constitute a complete set of differential 
invariants, and we will come to use these in place of the Wilczynski invariants from Section 4 
onwards. 
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2.2 Non-stretching variations 

Suppose that F : / x (— e,e) — t- is a smooth mapping such that, for fixed t, F(x,t) is 
a nondegenerate curve parametrized by x. Without loss of generality, we will assume that 
7(x) = F(x, 0) is parametrized by centroaffine arclength. Let X denote the variation of 7 in the 
t-direction, and expand 



dt 



F = a7 + 67' + c7". 



(We still use primes to denote derivatives with respect to x, although x is not necessarily an 
arclength parameter along the curves in the family for t ^ 0.) 

To compute the variation of the arclength differential |7, 7', 7"|^/^ dx, we introduce the no- 
tation 6 for variation in the f-direction along 7. Using the relation (2.2), we compute 

5Y = X' = (a' + poc)7 + {a + b'+ pic)l' + (6 + c')7", 

5-1" = X" = {a" + 2poc + p'qC + pob)l + {2a' + poc + b" + 2pic' + p[c + pib)l' 
+ {a + 2b' + pic + c")l" . 

Then 

(^|7, y, 7"| = |(^7, y, 7"| + |7, 57', 7"| + |7, 7', 6Y'\ = 3a + 36' + c" + 2pic. 
In particular, the variation X preserves the centroaffine arclength differential if and only if 

b' = -a-l{c" + 2pic), (2.4) 

i.e., 

X = al-(^j (a + y + fpic) dx^ 7' + c7". (2.5) 

We refer to vector fields of this form as non- stretching, since not only do such variations preserve 
the overall arclength of, say, a closed loop, but also no small portion of the curve is stretched or 
compressed. 



3 Hamiltonian curve flows 

3.1 Symplectic structure on starlike loops 

Generalizing Pinkall's setting [26] for planar star-shaped loops to the three-dimensional case, we 
introduce the infinite-dimensional space 

M = {7 : 5^ ^ : |7, Y, 7"| = 1}, 

as a subset of the vector space V = Map(S'^,M'^) of C°° maps from to M^. Assume that 
7 G M, i.e., 7 is a closed starlike curve parametrized by centroaffine arclength; then a vector 
field X = al + bl' + c7" is in the tangent space T^M if and only if X is of the form (2.5), where 
the coefficients a and c are 27r-periodic functions of x and satisfy the "zero mean" condition 

(a + |pic) dx = 0, 
ensuring that the coefficient of 7' in (2.5) is also periodic. 
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On V, define the skew-symmetric form 

uj^{X,Y) = (f \X,Y,Y\dx, X,YeT^V. (3.1) 

Note that uj is automatically closed (that is, doj = 0) since the integrand in (3.1) is a volume 
form on [1, 2]. 

Letting X = a7 + 67' + c7", Y = + b7' + c7", we compute 

uj^{X, Y)= i (ac - ac) dx. (3.2) 

Assuming that 1 and X,Y & T^M, then uj-y{X, Y) = for aU Y if and only if a = -|pic 

and c is constant. Thus, the restriction of a; to M is a degenerate closed 2-form (a pre-symplectic 
form), with kernel given by the subspace MZq + MZi of constant-coefficient linear combinations 
of the vector fields Zq = 7' and Zi = l" — |pi7 (corresponding to c = and c = 1 respectively). 

Note that this degeneracy is the result of restricting the 2-form (3.1) to the space of closed 
curves satisfying the arclength constraint. Degeneracy coming from constraints is common 
when defining symplectic structures on loop spaces [19] and phase spaces of nonlinear evolution 
equations [7]. 

Remark 3.1. The generators Zq and Z\ of the kernel of a; will turn out to be the seeds of the 
double hierarchy of curve flows discussed in Section 5. A similar situation is encountered for 
the LIE hierarchy [14], where the seed 7' spans the kernel of the natural pre-symplectic form on 
loops in Euclidean M^. 

One could attempt to remove the degeneracy by constructing a quotient of M with respect 
to group actions generated by flowing by Zq and Z\. Flow by Zq generates an action of the 
additive group M, simply by translation in x. The resulting quotient space M = M/M can be 
identified with the space of unparametrized starlike loops. Because of translation invariance, cj 
descends to a give a well-defined closed 2-form on M, with one-dimensional kernel MZi. On the 
other hand, we do not know of a natural geometric interpretation for the quotient of M by flow 
under on which w would become non-degenerate. 

However, can still be used to define a link between vector fields and functionals, and we 
will see that Z\ is linked in this way to the arclength functional. 



3.2 Examples 

Recall that the correspondence between vector fields Xu and (differentials of) Hamiltonians 
K € C°°(M) on a manifold M with symplectic form cj is defined by the relation 

d}i\X\ = LOy{X,XH), VX G TyM, (3.3) 

Xh being the Hamiltonian vector field corresponding to H. However, when w is degenerate 
the correspondence is no longer an isomorphism: for those functionals H for which there is 
a Hamiltonian vector field, Xh is only defined up to addition of elements in the kernel of oj. 

We will use (3.3) to compute Hamiltonian vector fields for a few interesting functionals; to do 
so, we will initially work in the ambient space V, and use the arclength-preserving condition (2.4) 
to rewrite the differential in a form suitable for applying (3.3). 

We first consider the arclength functional 

L(7) = /|7,7',7'f/^dx. (3.4) 
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on the space V. Given an arbitrary vector field X = al + hi' + cl" (not necessarily arclength 
preserving), the variation of the determinant in (3.4) along X is given by 

5|7,7',7"| = \X,l\l"\ + |7,X',7"| + |7,7',X"| = (3a + 36' + c" + 2pic)|7, Y, 7"|. 
Assume now that 7 G M, so that |7,7',7"| = 1 and 
dL[X] = i {a+ |pic) dx. 

We now seek a vector field Xl = al + hi' + cl" £ T^M such that dL[X] = oj^{X,Xl) = 
§^{ac — ac)dx. Using the non-stretching condition (2.4), we obtain the following Hamiltonian 
vector field 

Xl = Zi= 1" - Ivil 

(which is unique only up to adding a constant times Zq). In Section 4.4 we will see that the 
associated curve flow It = Xl leads to the Boussinesq equation for the curvatures ki, /c2- 
Next, we introduce the total curvature functional 

P(7) = (fpidx, 7 G M. (3.5) 



From 1"' = Pol + Pil' and (2.1), it follows that pi = |7,7"',7"|. Then the variation of pi along 
an arbitrary vector field X = al + hi' + cl" is given by 

5pi = |X,7"',7"| + \1,X'",1"\ + |7,7"',X"| 

= \X,pol +pil',l"\ + |7,X"',7"| + \l,pol + Pil',X"\ 

= Spia + 3c'po + 2cpo + 3a" + h'" + 4c"pi + Sc'p'i + cp'[ + bh'pi + hp'i + 2cp\. (3.6) 
Then, up to perfect derivatives, 

dP[X] = <T> 5pi dx = (p 3pia + 3c' po + 2cp'q + Ac'pi + 3cVi + cp'[ + Ah'pi + 2cp1 dx. 

J 'y J ^ 

Assuming X is an arclength-preserving vector field, we set h' = — a— |(c"+2pic) and compute 

® 5pi dx = (p —api + 3c' pq + 2cpo + |c"pi + 3c'p'i + cp'[ — \cp\ dx. 
J -^f J 

Integrating by parts, we arrive at 

dP{X) = <f -pia + {-p'q + y; - lpl)cdx. (3.7) 

Suppose that Xp = a7 + hi' + cl" is also arclength-preserving. Setting the right-hand side 
of (3.7) equal to a;^(X, Xp) = (p {ac — dc) dx, we get d = p'q — + |pf and c = —pi. Using 

equation (2.4) we compute b' = —{p'o — \pi + — g(— Pi — 2pf) = {p'l — po) , a perfect 
derivative. Thus, a Hamiltonian vector field corresponding to (3.5) is 

= (i {Pi - Pi) +Po)l + (Pi - Po)l' - Pil"- 

(This is only unique up to adding a linear combination of Zq and Zi.) Again, we will see that 
the associated curve flow It = Xp is also directly related, at the level of the curvatures, to one 
of the flows in the Boussinesq hierarchy. 
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4 Integrable centroaffine curve flows 

In this section we will examine the evolution of centroaffine curvatures induced by the curve 
flows defined in Section 3.2. We begin by computing the evolution of invariants under more 
general curve flows. 



4.1 Evolution of invariants 

First, we consider how the centroaffine invariants of a starlike curve evolve under a general 
non-stretching evolution equation 

7t = ro7 + riY + raY' (4.1) 

with ro = —r[ — ^(r2 + 2pir2). (Thus, we assume from now on that 7(x,i) is parametrized 
by arclength x at each time.) Of course, in order for (4.1) to represent a geometric evolution 
equation, ri and r2 should be functions of the invariants po,pi and their arclength derivatives. 

Proposition 4.1. The evolution equations induced by (4.1) for the Wilzcynski invariants are 

{Po)t = -r'i" + pir'l + 3por[ + p'^ri 

- I {4"' + Pi4') + ((3P0 - 2p;)r^)' + i {plr'2 + {p,p[ - p'/O-a) + p'^r^, (4.2) 
iPi)t = -2r'/' + 2pir[ + p[ri - r'^' + pir^ + (3y>o - Pa)^ + (Vo - PiO'^a- (4.3) 



Proof. The second equation (4.3) follows by substituting a = —r[ — ^(^2 + 2pir2), c = r2 in 
the last line of (3.6). Similarly, using po = |7"',7',7"|, we obtain (4.2) by computing 

{Po)t = |(7t)"',7',7"| + |po7+R7',(Y)',7"| + Po|7, 7', (Y)"|. ■ 

We note that these evolution equations previously appeared in [4]. 

From now on, we will take ki = pi and k2 = Pq — p'l as fundamental invariants; one reason 
for doing this is that the evolution equations for these invariants induced by (4.1) take the form 

(4.4) 



where V is the skew-adjoint matrix differential operator 

/ -2D^ + Dki + kiD -D^ + D^ki + 2Dk2 + k2D \ 

~ \ - kiD^ + 2k2D + Dk2 I {D^ + kiDki - kiD^ - D^ki) + [k2, D^] ) ' ^ ' ' 

D stands for the derivative with respect to x and [•,•] denotes the commutator on pairs of 
operators^. This operator V, which arises naturally when using ki, k2 instead of po; Pi^ will 
play a significant role in the integrable structure of the flows we study. 



4.2 Two integrable flows 

The vector field Xl induces a non-stretching evolution equation 

7t = 7" - Ikil. (4.6) 

^Note that expressions like Dki and Dk2 denote composition of D with multiphcation by ki and k2, respectively. 
The skew-adjointness of V is easy to check, given that D is skew-adjoint. 
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(This will be the first non-trivial curve evolution in the hierarchy discussed in Section 5, where 
the right-hand side is labeled as Zi.) By setting ri = 0, r2 = 1 in (4.4), we obtain the 
corresponding curvature evolution 



This PDE system for curvatures is Hamiltonian, since it can be written in the form 



where E denotes the vector-valued Euler operator 

on scalar functions / of fci, k2 and their higher x-derivatives kif^ ^ ki"^ . (One can check that 
the Poisson bracket defined using the Hamiltonian operator V on the appropriate function 
space - see Section 5.2 below - satisfies the usual requirements of skew-symmetry and the 
Jacobi identity.) 

Moreover, (4.7) can also be written in Hamiltonian form as 

^j)^ = QEp3, (4.9) 
for a different Hamiltonian operator and density 



Q=(° o), Pz:=W:i' + k2k!^^kl^\k\. (4.10) 



(The notation pj, is explained below.) Since the curvature evolution can be written in Hamil- 
tonian form in two ways (4.7) and (4.9), the integrals J k^ dx and J psdx are conserved by the 
flow (for appropriate boundary conditions). 

Remark 4.2. In fact, the curvature evolution here is a bi- Hamiltonian system, because V and Q 
are a Hamiltonian pair, i.e., their linear combinations form a pencil of Hamiltonian operators, 
and a pencil of compatible Poisson structures. This assertion can be verified mechanically (see, 
e.g.. Section 7.1 in [22] for details), but it also follows from the fact that, at least in the periodic 
case, the Poisson structures are reductions of a well-known compatible pencil of Poisson brackets 
on the space of loops in s[(3). (Indeed, when 7 is periodic -or more generally has monodromy - 
the matrix K in (2.3) provides a lift into this loop space.) The proof of the reduction of these 
brackets can be found in [8], where V is linked to the Adler-Gel'fand-Dikii bracket for s[(3) 
and Q is associated to its companion. The brackets were later linked to curve evolutions and 
differential invariants in [15], where more details are available. 

The (negative of the) Hamiltonian vector field of Section 3.2 induces the non-stretching 
evolution 

It = kil" + k2l' + rol, (4.11) 

where 

ro = -{k'^ + \[k'i + 2kl)). (4.12) 
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(The right-hand side of (4.11) is labeled as Z2 in the hierarchy discussed in Section 5.) We 
similarly obtain the curvature evolution equations induced by this flow by setting r\ = A;2, 
r2 = h\ in (4.4). We remark that the resulting system is also bi-Hamiltonian, since it can be 
written as 



PEp2 = QE/j4 (4.13) 



for 



P2 = hk2, PA = \{k'lf + k'l{k'2 - kj) - ki{k[f + (4)2 - klk'2 + lk\ + 2kikl 



Thus, j p2 dx and J p^dx are conserved integrals for (4.11). (Because Xp corresponds symplec- 
tically to the Hamiltonian J ki dx, it is automatic that this integral is also conserved.) 

Remark 4.3. The arclength normalization (2.1) is preserved by the simultaneous rescaling 
X ^ Xx, 7 I—)- A^^7. Under this rescaling, ki and ^2 scale by multiples and respectively. 
Thus, we may assign scaling weights 2 and 3 respectively to these curvatures, and each x- 
derivative taken increases weight by one. 

It will turn out (see Section 5.1 below) that the conserved densities for evolution equa- 
tions (4.6) and (4.11) are all of homogeneous weight, with one density for each positive weight 
not congruent to 1 modulo 3. We will number the densities in order of increasing weight, letting 
Po = ki, pi = k2 and so on; thus, the density in (4.9) is denoted by /?3, since its weight falls 
between those of p2 and p^. 

The curve flows (4.6) and (4.11) turn out to share the same conservation laws; for example, 
/ ki dx is conserved by (4.6) because (4.7) implies that 

{k^)t = D{k[ + 2k2). 

Similarly, (4.11) conserves / k2dx because (4.13) implies that 

{k2)t = D (lk[^^ + k2 - 2kik'l - {k[f - 2kik'2 + ^kl + 2k^ . 

In Section 5 we will show that these flows share an infinite sequence of conservation laws. 



4.3 Lax representation 

In this subsection we use geometric considerations to derive Lax pairs for curvature evolution 
equations induced by (4.6) and (4.11). 

In [3] we found that the components of the solution 7(2;, t) of Pinkall's flow satisfied the 
scalar Lax pair for the KdV equation. In the same spirit, we seek a system of the form 

Cy = 0, yt = My, (4.14) 

satisfied by each component of 7, where C and M are differential operators in x with coefficients 
involving ki, k2- Using (2.2), we see that every component of 7 satisfies the scalar ODE y'" = 
(f^iy)' + k2y, and so we will let 

C:=D^- Dki - k2 

and seek operators Aii for (4.6) and M.2 for (4.11). 

In the case of (4.6), the components of 7 also satisfy yt = y" — ^kiy, so we choose 

Ml := - ffei. 
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One can then verify that (4.7) implies that 

Ct = [MuC\. (4.15) 

In the case of (4.11), the components of 7 satisfy yt = kiy" + k2y' + roy, with ro as given 
by (4.12). So, we might set M.2 = kiD'^ + /c2-D + tq. However, (4.14) would also be satisfied 
if we modify A^2 by adding J\fC, where M is an arbitrary differential operator. In fact, the 
system (4.13) actually implies that Ct = [7^2) -C] for 

M2 := {kiD^ + k2D + tq) - 3DC. 

Writing these systems in Lax form (4.14) enables us to interpolate a spectral parameter into the 
linear equations satisfied by the components. Thus, consider solutions of the compatible system 

Cy = Ay, yt = Mjy, (4.16) 

where j = 1 or j = 2. Of course, the components of the evolving curve satisfy (4.16) only when 
A = 0. When A 7^ 0, we can construct solutions of the curve flow using solutions of (4.16): 

Proposition 4.4. Let ki, k2 satisfy the evolution equation (4.7) for j = \ or (4.13) for j = 2. 
For fixed A S M, let yi, 2/2, 2/3 be linearly independent solutions of (4.16), with Wronskian W. 
Then W is constant in x and t, and 7 = W~^/^{yi^y2-,y^)^ is arclength-parametrized at each 
time t, with centroaffine invariants ki and k2 = k2 + A. Furthermore, 7 satisfies the evolution 
equation 



It 



7" - Ikil, j = 1, 

kil" + {k2 - 4A)r + ro7, j = 2. 



Proof. If we let y = (yi, 7/2,2/3)"'" and form the matrix F = (y, y', y"), then F satisfies differential 
equations of the form 




F-'F^ =110 ki I , F-^Ft = N„ 



where both right-hand side matrices have trace zero. For example, when j = 1 one can directly 
calculate, by differentiating yt = M.iy, that 



f-jki k2 + lk[ + X k'^ + \k'i 

\ki k2 + lk[ + A 

V 1 ifci 



Thus, the Wronskian W is constant in x and t. 

Because 7"' = (kil)' + (/c2 + A)7, the centroaffine invariants of 7 are ki and ^2- It is 
straightforward to compute 7f in the j = I case, using yt = A4iy. In the j = 2 case, we compute 

yt = M2 y = kiy" + {k2 - \)y' + r^y - 3DCy = kiy" + {k2 - 4A)y' + roy. ■ 



4.4 Connection with Boussinesq equations 

In [4] Chou and Qu note that, under the centroaffine curve flow (4.6), the curvatures ki, k2 
satisfy a two-component system of evolution equations that is equivalent to the Boussinesq 
equation. This suggests that the other integrable flow (4.11) under discussion may be related 
to the Boussinesq hierarchy. 
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Dickson et al. [5] write the (first) Boussinesq equation as a system 

{qo)t + yi +yiqi=0, {qi)t-2q'Q = 0. (4.17) 
They embed this in a hierarchy of integrable equations, each of which is written in Lax form as 

Lt = [P^, L], L:=D^ + qiD + y[ + qq, (4.18) 

where Pm is a differential operator of order m ^ mod 3, with coefficients depending on goj Qi 
and their x-derivatives. Note that Pm must be chosen so that [Pm, L] has order one. For example, 
while Pi = D yields the trivial evolution {qi)t = q[, {qo)t = Qq, setting P2 = + |qi gives the 
Boussinesq equation (4.17). 

Given the resemblance between (4.15) and (4.18), it is tempting to find substitutions to 
connect the Boussinesq equation with (4.7). In fact, we can make L and C coincide by setting 

fci = -gi, k2 = yi-qo. (4.19) 

With this substitution, Mi coincides with P2, so it follows that (4.7) and (4.17) are equivalent. 

In [5] it is shown how the coefficients of the operators Pm can be obtained solving a recursive 
system of differential equations, and thus these depend on a number of constants of integration. 
For example, the expression for P4 is 

Pa = [fiD^ + {91 - y[)D + (i/f -g'i + |gi/i)] + 

+ [foD^ + {90 - \QD + ilf^ -g'^ + |gi/o)] L + k^^ + ^iL, 

where /o = 0, go = 1) /i = g9i + ci, gi = \qQ + di, and /c4^0i ^4,ii ci, di are arbitrary constants. 
For convenience, we will set all these arbitrary constant to zero, so that 

P^ = D^ + jqiD' + liq'i + qo)D + + jq'^ + jql 

Again, if we use the substitutions (4.19), we find that the operator A^2 coincides with — 3P4. 
Thus, (4.13) is equivalent to the second nontrivial flow in the Boussinesq hierarchy, provided we 
also rescale time by t — t- — 3t. 



5 Hierarchies 

In [22] the Boussinesq hierarchy is discussed as an example of a bi-Hamiltonian system, in 
which two sequences of commuting flows (and conservation laws) are generated by applying 
recursion operators. Thus, given the equivalences established in Section 4.4, it is not surprising 
that the Poisson operators defined in Section 4 can be combined to give a recursion operator 
that generates a double hierarchy of commuting evolution equations for ki, k2- In fact, we will 
show that our recursion operator is equivalent to the Boussinesq recursion operator as given in 
Example 7.28 of [22]. The new information we add is that each of these evolution equations is 
induced by a centroaffine geometric evolution equation for curves, which is itself Hamiltonian 
relative to the pre-symplectic structure defined in Section 3.1 (see Theorem 5.4 below). 

5.1 Recursion operators 

We define a sequence of evolution equations for ki , /c2 



(5.1) 
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via the recursion 

Fj+2 = VQ-^Fj, (5.2) 
with initial data given by 

p — ( k'{ -\- 



Vs^^l^l ^1 ) ^2 

(Note that Fi is the right-hand side of (4.7), while for j = (5.1) gives a simple transport 
equation for ki, k2, corresponding to flow in the direction of the tangent vector 7'.) 

In order to assert that the Fj defined by (5.2) are local functions of ki, /c2 and their deriva- 
tives - i.e., in calculating each Fj, the operator is only applied to exact x-derivatives of local 
functions - we cite well-known results on the Boussinesq hierarchy. For example, the version of 
the first Boussinesq equation used by Olver [22] is 

Ur = v' , Vr = \u"' + (5-3) 

where r is the time variable. If one considers linear transformations on the variables, it is 
necessary to use some imaginary coefficients to make our version (4.7) of the first Boussinesq 
equation for ki, k2 equivalent to (5.3): 

X = X, T = it, ki = —2u, k2 = u' — iv. (5-4) 

Proposition 5.1. Under the above change of variables, the recursion operator VQ"^ is equiv- 
alent to the Boussinesq recursion operator in [22] . 

Proof. The transformation between ki, k2 and u, v can be written as 
'ki\ ^ ^ f-2 



k2j \vr ^ \D 

Thus, if d/dt{u,v)^ = F[u,v] is an evolution equation for u, v, the right-hand side of the 
corresponding evolution for fci, A;2 is ^ o F. Thus, our recursion operator VQ~^ for flows on 
the ki, k2 variables corresponds to a recursion operator 

g-^VQ-^g (5.5) 

on flows in the u, v variables. In fact, when one calculates (5.5) and substitutes for ki, k2 in 
terms of u, v using (5.4), the result is exactly — i times the Boussinesq recursion operator given 
in [22]. ■ 

Since in [27] (see Section 5.4 in that paper) it is proven that the Boussinesq recursion operator 
from [22] always produces local flows when applied to the 'seed' evolution equations (i.e., the 
tangent flow and first Boussinesq), it follows that the same is true for our recursion operator. 

Remark 5.2. Once one checks that the evolution equations (5.1) for j = and j = 1 commute, it 
is automatic from the bi-Hamiltonian structure that all evolution equations in the sequence (5.1) 
commute in pairs (see, e.g.. Theorem 7.24 in [22]). 

It is easy to check that the 'seeds' Fq, Fi for the recursion are related to the initial conserved 
densities by 



Fo = VEpo = QEp2, Fi = VEpi = QEp^. 



(5.6) 
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Table 1. 



Non-stretching vector field 


Conserved density 


Zo = l' 


Po = ki 


Zi = 7" - Ikil 


Pi = k2 


Z2 = kil" + k2l' + ... 


P2 = kik2 


Zs = {k[ + 2^2)7" + {\kl - Ik'i -k'2)r + ... 


P^ = \{K? + k'^k2 + lkl + kl 


Z4 = {-k'C - 2k'^ + 2kik[ + 4kik2p" + (f/ti^^ + k'^' 
-2kik'l - {k[f - 2kik'2 + ^kf + 2A;i)7' + . . . 


Pi = Wif + k'i{k'2-kl)-k^{k[f 

"^(^2} ^1^2 9^1 '2i}\i\l\i<2 



(The coefficient of 7 in some vector fields is omitted for reasons of space, but can be determined from 
the non-stretching condition.) 

• The 7' and 7" coefficients of Zj match the components of Epj. 

• Densities satisfy the recursion relation = Q^^VEpj 

• 7f = Zj induces curvature evolution | ^ ) = "PEp, = QEp, 1 2. 

\k2Jt 

• For j > 2, Zj is a Hamiltonian vector field for — J Pj-2 dx. 



(The second set of equations was derived in Section 4.2.) While VQ~^ is the recursion operator 
for commuting flows, it is evident from (5.6) that Q^^V should be the recursion operator for 
conservation law characteristics (i.e., the result of applying the Euler operator E to a density). 
In fact, we may define an infinite sequence of conserved densities by 

Epj+2 = Q-^VEpj, j > 0. (5.7) 

The first few densities calculated using this recursion appear in Table 1. 

We now use these densities to define a sequence of flows for centroaffine curves, and relate 
each of them to a curvature evolution equation in the sequence (5.1). Namely, if / is any local 
function of ki, k2 and their derivatives, we define 

X/ := (E/)i7' + (E/)27" + ro7, (5.8) 

where the subscripts indicate the components given by (4.8) and ro is determined by the non- 
stretching condition. Then for the sequence of densities defined recursively by (5.7) we define 
the vector fields 

Zj := (5.9) 

and the corresponding sequence of curve flows 

It = Zj. (5.10) 

Proposition 5.3. For each j > the curve flow (5.10) induces the curvature evolution 
g{kiM? = Fj. 

Proof. From (5.6) and the recursion relations, it follows by induction that 
Fj=VEpj, j>0. 



Then the result follows immediately from (4.4). 
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5.2 Hamiltonian structure at the curve level 

We now consider the question of how the Hamiltonian operator V is related to the Hamilto- 
nian structure defined at the curve level in Section 3.1. Recall from Section 3.2 that X is 
a Hamiltonian vector field associated to the functional H if 

dH[Y] =uj^{Y,X) 

for any non-stretching vector field Y . 

Theorem 5.4. Let H{1) = (p pdx and assume that 

Ep=Q-^rEp, (5.11) 

i.e., p is next after p in the sequence of densities generated by the recursion operator Q^^V. 
Then —X^ (as defined by (5.8)) is Hamiltonian for H . 

Proof. Based on the definition (3.1) of cj, we need to show that dH[Y] equals 

^\Y,y,-XP\dx = j\XP,-i\Y\dx yYGT^M. 

U X = al + 67' + c7" and F = a7 + 67' + c7", then from (3.2), 

(j) \X, 7', Y\dx= i> {ac - ac) dx. 
J ^ J ^ 

However, using (2.4) to eliminate a and a, we obtain 

j \X, 7', Y\dx = j (^-b'c + cb' + \{cc - c"c)) dx = - j{b'c + c'6) dx, 
where the last equation follows by integration by parts. Thus, 

\XP,l',Y\dx = - / r ) -QEpdx = - / r ) -V^pdx, 



7 ^7 \v J7 \C/ 

using (5.11) in the last step. Then, because V is skew-adjoint, 

j \XP,l\Y\dx = jEp-v(^^ dx. 
On the other hand, using the properties of the Euler operator we have 

dH[Y] = (f{Ep)i5Yh + (Ep)25yfc2dx, 
where 5y denotes the first variation in the direction of Y. Now using (4.4) we have 

This concludes the proof. 

The following corollaries are immediate consequences of the theorem. 
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Corollary 5.5. Define the Poisson bracket 

{H, G} = j) Eh-VEgdx, 

where G{ki,k2) = § gdx and H{ki, k2) = j>hdx, and g, h are functions of periodic ki, k2 and 
their derivatives. Then dH[X^] = {H,G}. 

Corollary 5.6. The vector fields Zj defined by (5.9) are Hamiltonian for j > 2. 

Corollary 5.7. A closed curve 7 is critical for the functional H{1) = <f^ pj dx with respect to 
non- stretching variations if and only ifl is stationary for a constant- coefficient linear combina- 
tion Zj^2 + CqZq + CiZi. 

Proof. By Theorem 5.4, —Zj^2 is Hamiltonian for H. Thus, a curve 7 is i?-critical if and only 
if uj{Y, Zjj^2) = for all Y G T^M. This condition is satisfied if and only if Zj j^2 is in the kernel 
of uj^, i.e., along 7 it is equal to a constant-coefficient linear combination of vectors Zq and Zi. 
Equivalently, + co^q + ciZi = along 7 for some constants cq, ci, expressing the property 
that 7 is stationary for a linear combination of these vector fields. ■ 



5.3 Projective properties 



As stated in Remark 2.2, a centroaffine curve 7 projects to give a conic in MP if and only if the 
Wilczynski invariants satisfy pq — ^^p'l = 0. (The corresponding condition in terms of fci, /c2 is 
^2 + 2^1 = 0.) In this subsection we will investigate flows in the hierarchy having the property 
that, if 7 projects to a conic at time zero, then it continues to have a conical projection at 
subsequent times. We will show later that the equation of the conic in homogeneous coordinates 
is fixed in time. We will also discuss flows that preserve a parametrization that is proportional 
to projective arclength; in that case, the corresponding condition in terms of curvatures is that 
^2 + is a nonzero constant along the curve. 

These investigations are much easier if, instead of k2, we use an invariant that vanishes 
precisely when the condition we are investigating holds. Accordingly, we fix a constant C, and 
define an alternative pair of invariants 



u 



ki, v = k2 + Ik'i - C. (5.12) 



Thus, the curve is a conic if v = when C = 0, and the curve has a constant-speed parametriza- 
tion (relative to projective arclength) if u = when C 7^ 0. 

We will convert the evolution equations in the hierarchy (at the level of the invariants) to 
these variables. Suppose that a curve flow causes the invariants to evolve by 



ll^^ = F[k,,k2], 



where -F is a vector- valued function of ki, k2 and their derivatives. Then the corresponding 
evolution equation for the alternative invariants is 



\D 1 



Similarly, if IZ is the recursion operator generating the hierarchy of evolution equations 
for fci, k2, then the recursion operator for the corresponding flows on n, v differs by a gauge 
transformation: 



'R, = QnQ-^, where Q'^ 



1 0^ 

-\D ir 
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(It is understood that, in TZ on the right-hand side, ki, k2 are substituted for in terms of u, v.) 
Specifically, using TZ = VQ^^ as defined by (4.5), (4.10), we compute 

n = TZo + FoD-^ {0 l)+FiD-^{l O) , 

where 

^(v + C) 2[u-D^)\ ~ _ fu'\ ~ _f 2v' 



and Af is the scalar differential operator gZ)^ — ^uD^ — \u' D + — \u" . One can check that 
the vectors Fq, Fi are the time derivatives of n, v, corresponding to the 'seeds' Zq and Zi for 
the hierarchy of curve flows. 

By applying IZ to Fq, Fi, one can generate the right-hand sides of the evolution equations in 
the hierarchy in terms of u and v. Letting Fj denote these vectors, we compute (for example) that 



F2 = IZFq 



-2v"' + A(uv)' + -iCv! \ 
1^(5) _ ^^/// _ 2v!v!' + \v?v! + (4r; + ?,C)v') ' 

+ \(v?v! - uv!") - fu'u" + (10?; + QC)v' \ 
i^;(5) _ (5^ + ic)u"' + l{{u^v)' - uv'" - u"v') - \v!v" + 2CWJ ' 



Here, when applying to differential polynomials in -u, u , the constant of integration is taken 
to be zero. 

Notice in particular that if u = 0, then the bottom component of -F3 — "iCFx vanishes. Thus, 
the flow Z3 — 3CZi preserves the condition that v is identically zero. In fact, we can calculate 
two infinite sequences of evolution equations for n, v that preserve this condition; the right-hand 
sides of these are 



k 



^Q){-3CyF2ik-j), fceven. 



(5-13) 

^(5)(-3C)^F2(fc-i)+i, fcodd. 



While it is routine to verify that any individual curvature evolution equation in these sequences 
preserves = 0, it is easier to observe that the members of these sequences satisfy the recursion 
relation Gk+2 = (^ ~ "^CYGk- Then the fact that they all preserve v = is a consequence of 
the following: 

Proposition 5.8. If a curvature evolution (ut,vt)'^ = Gk[u,v] in this sequence preserves v = 0, 
then so does the evolution {ut^vtf^ = Gk+2[u,v]. 

Proof. We assume that Gk = {F>ii, Dl2)^ for local functions £1, £2 of u, v and their derivatives. 
(This form for Gk is necessary if we are able to apply operator IZ to it and produce local 
functions.) 

Within the ring of polynomials in n, v and their derivatives, let V denote the ideal generated 
by v, v' , v" , etc. By hypothesis, DI2 G and the same is true for ^2- 
We compute 
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Thus, the bottom component of {1Z — 3C)'^Gk is given by 

{M + Fi2D-^) {iiFu + I2F01 + 3vDii + 2{u - D^)Di2) 
+ {3v + Fo2D-^) {£iFi2 + £2Fo2 + NDh + 'ivD£2) ■ 

(Here, Fji and Fj2 denote the top and bottom entries in the vector Fj.) The second factor in 
the top hne is the top entry of {IZ — 3C)Gk- This polynomial lies in V, and the same is true if 
we apply A/" or to it. On the other hand, because F02 = v' £ V, the coefficient in front on 
the second line also vanishes when v = 0. ■ 



In the special case when C 
preserve conicity: 



0, we see that the following curve flows (as defined by (5.9)) 



Z3, Z^, Zj, Zg, Zii, z 



12, 



(5.15) 



However, when C 7^ 0, some care needs to be taken in matching the evolution equations for u, v 
that preserve v = with the corresponding linear combinations of the curve flows Zj. 

In the proof of Proposition 5.8, we used the fact that if an exact derivative D£ lies in V, then 
by choosing the constant of integration equal to zero, the antiderivative i also lies in V. However, 
if we express D£ in terms of ki, k2 instead of u, v, and then take an antiderivative, a particular 
constant of integration must be chosen in order to belong in V. Thus, when we compute the A;th 
evolution equation for u, v by using the recursion operator IZ (which involves applying D^^), 
then convert this to an evolution equation for ki, k2, and finally try to match it with a curve 
flow in the hierarchy (5.10), we get a linear combination of with lower-order flows of the same 
parity. For example, if we substitute (5.12) into F2, and then apply the operator Q~^, we get 



g-^F2[ki,k2+lk[-c] 



2ui 
,3^1 



F2 — CFt) 



-k'{" 



2k'^' + 2kik'l+4:{kik2y +2{k[f 
2kik'l'-4k[k'l+4k2k'2 



-Ck[ 



2ikik'^y + l{ki)%- 



' Ckrf 



Thus, F2 is induced by the curve flow Z2 — CZQ] similarly, F^ is induced by Z3 — 2CZi, and so on. 

Similarly, when we apply the recursion operator iJZ — 3C)^ to generate higher-order flows 
that preserve u = 0, these constants of integration accumulate and change the relatively nice 
pattern of the coefficients exhibited by (5.13). Here is what we get when we compute the first 
few curve evolutions corresponding to the evolutions Gk'- 



u, V evolution 


centroaffine curve flow 


Go 
Gi 
G2 
G3 
G4 


Zo 

Z3 - 5GZi 

Z4 - 7CZ2 + UC^Zo 

Zj - IICZ5 + AAC^Zs - ^C^Zi 

Zg - ISGZe + QbC^Zi - ^C^Za + ^G^Zq 



5.4 Conical evolutions and the Kaup— Kuperschmidt hierarchy 

In this section we examine special properties of the conicity-preserving flows (5.15) (hence, from 
now on we are assuming C = 0). These properties will enable us to connect our hierarchy of 
centroaffine curve flows in with the Kaup-Kuperschmidt hierarchy and with curve flows 
in centroaffine M^. We begin with the observation that, when restricted to conical curves, the 
coefficient of 7" vanishes for as many of the vector fields in (5.15) as one cares to check. In other 
words, this coefficient belongs to V, the ideal within the ring of differential polynomials in ki, k2 
generated by k2 + ^k[ and its derivatives. In fact, this is true in general, as shown in the following: 
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Proposition 5.9. If j = or j = —1 modulo A, then the bottom component ofEpj belongs in V; 
hence, the 7" coefficient of Zj vanishes on conical curves. 

Proof. The statement can be verified directly for j = and j = 3. For higher values, we use the 
recursion relation between the characteristics, which implies that Epj^4 = {Q~^T')'^Epj. From 
Proposition 5.8 we know that the bottom component of Fj = QVEpj lies in V. By inserting 
powers of Q and into the recursion relation, we get 

Epj+4 = Q-^vQ-^g-^gvEpj = Q-^izg-^Fj. 

As in the proof of Proposition 5.8, we can assume that Fj = {Dii,Di2f where £2 G V. Taking 
C = in equation (5.14), we see that the top entry of TZG'^Fj = Q'^TZFj is 

hFn + hFoi + ^vDli + 2(n - D'^)D^2, (5.16) 

which clearly is in "V. (Note that v = k2 + \k'^ here.) Noting the form of Q, we see that the 
bottom entry of Epj^^ is applied to (5.16), so it must also belong to V. ■ 

Next, we make a connection with curve flows in centroaffine M? to show that, for the flows 
in (5.15), the cone that the curve lies on is preserved by the time evolution. 

Proposition 5.10. Ifl{x,t) evolves by any of the vector fields in (5.15), and 7(x,0) lies on 
a cone through the origin in M^, then l{x,t) lies on the same cone at later times. 

Proof. Using the action of SL(3) we can, without loss of generality, assume that the equation 
of the cone is yiy^ — (2/2)^ = 0. We fix a map V from onto this cone: 

Of course, when we projectivize on each end this gives the Veronese embedding of MP^ as 
a quadric in MP^. The scale factor of 2~^/'^ is chosen so that if X{x) is a parametrized curve in 
satisfying the centroaffine normalization |X, X'| = 1, then F = VoX satisfies the normalization 
|r,r',F"| = 1. Moreover, if p{x) is the curvature of X, then the invariants of F are ki = —4p 
and k2 = 2p'. Finally, if X evolves by the non-stretching flow 

Xt = rX' - \rX, (5.17) 

then F(x, t) = V o X(x^ t) satisfies F^ = rF' — r'F. 

By Proposition 5.9 all the flows in (5.15), when restricted to conical curves, are of this form, 
for some choice of differential polynomial r in ki. For any initial data 7(x,0), we can define 
a curve X{x, 0) in such that 7(x, 0) = \/oX{x, 0) and make X{x, t) evolve by (5.17). Because 
F(x,t) = V o X{x,t) satisfies the same initial value problem, then l(x,t) = T{x,t) at all times, 
and 1{x, t) lies on the cone defined by yiys — (1/2)^ = at all times. ■ 

In [4] Chou and Qu discovered a non-stretching flow for curves in centroaffine M'^ which 
preserves the conicity condition k2 + ^k[ =0 and which causes the curvature ki to evolve by 
the Kaup-Kuperschmidt equation: 

///// r /// 25 / // I r 2 / 

ut = u — buu — u + ou u 



(see Case 3 in Section 3 of their paper, taking A = 0). In fact, up to a multiplicative factor 
of 1/3, Chou and Qu's flow is the same as the restriction of Z3 to conical curves. 
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Not only does flow Z3 give a geometric realization of the Kaup-Kuperschmidt equation, but 
the entire sequence (5.15) of flows realizes the Kaup-Kuperschmidt hierarchy, when restricted to 
conical curves. To see this, note that the square of the recursion operator TZ relates the evolution 
of ki = u under Zj to its evolution under (Here, we use the notation of Section 5.3 but 

with C = and v = because of conicity.) The resulting recursion operator is 

-iD^^ + 2uD^ + 6n'Z)3 + [fu"-3u^) + (f n'"-10W) D + fu"" - fuu" 
- fiu'f + lu' + u'D-' o ilu^-lu") + I {u""'-5uu"'-fu'u"+5u\') D'K 

This agrees with the known recursion operator for symmetries of the Kaup-Kuperschmidt hierar- 
chy. (See, e.g., Example 2.20 in [29], where the operator differs by changing u to —u and rescaling 
time by a factor of 1/3.) Using this, one can check that the curvature flows induced by (5.15) 
for conical curves coincide with the commuting flows of the Kaup-Kuperschmidt hierarchy. 

Remark 5.11. The curve flow discovered by Chou and Qu is in fact also defined for centroaffine 
curves parametrized proportional to projective arclength (i.e., those for which k2 + \k'i is a con- 
stant), and nevertheless still induces Kaup-Kuperschmidt evolution for ki. Recently, Musso [20] 
has extended this to a hierarchy of flows for arclength-parametrized curves in MP'^ which induces 
the Kaup-Kuperschmidt hierarchy for curvature evolution. We suspect that these flows coincide 
with the restrictions of the flows studied in Section 5.3 (for C 7^ 0) to the centroaffine lifts of 
such curves in MP^. 

Remark 5.12. Schwartz and Tabachnikov [28] showed that certain maps defined on the space 
of convex polygons preserve the subset of polygons that are inscribed (or circumscribed) on 
a conic: that is, if the vertices of the polygon (or those of its projective dual) lie on a conic, then 
the same is true for its image under the map. The building blocks for these maps are elementary 
maps Tr that associate to a given polygon another polygon obtained from the intersections of 
diagonals joining each vertex to the vertex located r positions to the left or right of it. In fact, 
the maps preserving conicity are particular combinations of Tr for certain values of r. 

The map corresponding to r = 2 is called the pentagram map and it is known to be a dis- 
cretization (in both time and space) of the Boussinesq equation [24]. It is natural to wonder if 
the maps in [28] are somehow associated to flows in the Boussinesq hierarchy. We are currently 
investigating this. 
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